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Abstract 

The interpretation of virtual gluons as ghosts in the non-Unear ghionic structure of QCD permits 
the formulation and realization of a manifestly gauge-invariant and Lorentz covariant theory of in- 
teracting quarks/anti-quarks, for all values of coupling. The simplest example of quark/anti-quark 
scattering in a high-energy, quenched, eikonal model at large coupling is shown to be expressible as 
a set of finite, local integrals which may be evaluated numerically; and before evaluation, it is clear 
that the result will be dependent only on, and is damped by increasing momentum transfer, while 
displaying physically-reasonable color dependence in a manner underlying the MIT Bag Model 
and an effective, asymptotic freedom. Similar but more complicated integrals will result from all 
possible gluonic-radiative corrections to this simplest eikonal model. Our results are compatible 
with an earlier, field-strength analysis of Reinhardt et al. 
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I. INTRODUCTION 



There has long been a strong-couphng framework in Abehan QFT, whose lowest-order 
approximation is the Eikonal Model; and, with due attention to color indices and their dis- 
ruptive effects on the coherence of Abelian eikonalization, these techniques can be extended 

n 

to QCD [1]. In this spirit, we would like to call attention to a novel, manifestly gauge- 
invariant (MGI) method of calculating the sum of all virtual-gluon exchange graphs in 
QCD, including — and, in fact, made possible by — cubic and quartic gluon interactions. We 
illustrate this technique by its application to quark-quark (QQ) or quark-anti-quark(QQ) 
scattering in an eikonal- style, quenched approximation; in effect, we concentrate on the sum- 
mation of all possible gluon exchanges with color coupling constant g treated as an averaged, 
or constant quantity, neglecting its renormalization, along with quark mass and propagator 
renormalizations. 

By 'eikonal model' we mean one specific restriction: that all virtual-gluon 4-momenta 
emitted or absorbed by the scattering quarks are small compared with the incident and final 
4-momenta of the quarks in their center of mass (CM). Corrections to such an eikonal model 

n 

were defined years ago [2|], and in principle may be adjoined to the present discussion; but 
that is outside the present analysis. In Abelian Physics, this assumption leads to coherent 
scalar or Neutral Vector Meson (NVM) exchanges; in QCD specific color fiuctuations are 
introduced which can destroy such coherence. These color techniques were first introduced 
as an intelligent, quasi-Abelian (QA) approximation to a theory of simple non-Abelian 
exchanges; but, in the present paper, with its emphasis on MGI and a concurrent summation 
over all cubic and quartic gluon interactions, such approximations become exact. 

We treat the quarks as effectively asymptotic particles, since it takes only two or three 
in combination to make an asymptotic hadron. One can continue to retain gluons in the 
formalism, and if needed, introduce Faddeev-Popov ghosts [4] to insure that, when the such 
gauge-dependent gluon propagators are renormalized, expected properties are maintained. 
But in this paper such considerations are suppressed, for we are concerned only with the 
enumeration and summation of all virtual gluon exchanges; and it is for these exchanges 
that the MGI properties hold. It should also be mentioned, and will be noted below at an 
appropriate point, that all possible gluon-exchange corrections to the relatively simple forms 
presented below will possess the MGI property. 
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Perhaps the most frequently-used way of introducing gauge invariance in QCD is by the 
use of the functional integral (Fl) 



Z\j,7],r]] = U / d[A]8[T[A)]d.ei[5Tl5uj\-e^Y> 

J r]-G,[A]-r] + L[A]+t J j-A 



(1) 



• exp 

where 

G,[A] = [m + ^-{d-igA-\)r\ 
L[A] = Trln[l -ic/7- A- ASc], 

and where j^, r^^, and f/^ are gluon and quark sources, respectively, the delta-functional of 
J-'[A] defines the particular gauge adopted, and the det[6J^/6u!°'] guarantees the color-gauge 
invariance of the FI when a change of gauge is made by the variation of a relevant function 
uj^{x). A/" is a normalization constant which is chosen such that ^[0,0,0] = 1, and the FIs 
over quark coordinates have already been performed. The enumeration of gluonic degrees 
of freedom in this formalism is muted, but perturbative expansions of ([1]) are equivalent to 
those obtained immediately below. 

There is another, independent method of arriving at the equivalent of ([1]) in which one 
starts from Schwinger's Action principle [5], where the enumeration of proper degrees of 
freedom is paramount, while gauge invariance takes a secondary and circuitous path. There, 
as in QED, one immediately finds that the equal-time-commutation-relations (ETCRs) of the 
gluon field operators lead to proper quantization in the Coulomb gauge; but because of the 
micro-causality of the fields, [A^, A^]\^^_^^ = 0, the more complicated, canonically-conjugate 
field momentum operator ir^-lx) may be replaced by OqA^ for purposes of calculating relevant 
propagators in a variety of gauges. 

Schwinger's formalism is the one that we shall initially adopt, beginning with the choice of 
a relativistic gauge {e.g., one of standard gauges used in QED, or an axial gauge) for the free 
gluon generating functional (GF), expressing the full GF as a well-defined Action operator 
acting upon the free GF for gluons and quarks, and then employing a convenient rearrange- 
ment of the functional operations in terms of an equivalent but conceptually-simpler linkage 
operator. For specific processes, in that selected gauge, with the aid of Halpern and Frad- 
kin/eikonal representations, one sums over ALL virtual, gluonic fluctuations, including those 
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due to cubic and quartic gluon interactions; and one then trusts that subsequent events pro- 
vide the necessary gauge invariance, at least for all physical processes, and renormalizations, 
as in QED 

In this presentation, we begin as above; but before all gluonic fluctuations are performed, 
we observe that, in QCD, unlike QED, there is one special way of insuring MGI. This 
simple step corresponds to treating virtual gluons as ghost gluons, with the result that all 
of the initial, gauge-dependent gluon propagators cancel away; this is gauge invariance with 
a vengeance! With one small exception — which will be discussed and justified in Section 
Hn — the form of this result has previously been found in field-strength analysis, e.g., that of 
Reinhardt et al. [6], with the difference that our result is gauge independent, while that of 
Ref. 3 allows the choice of an arbitrary gauge; this difference is discussed following Eq. (12 9 p 
of the present paper. It may also be of interest to note that the analysis of Ref. 6[ is in 
Euclidean space, while ours is directly in Minkowski space. And in the context of an eikonal 
model, the 'local' simplifications obtained are sufficient to reduce all the FIs (used in the 
Fradkin representation of Gc[v4] and L[y4]) to ordinary integrals, susceptible to numerical 
integration. 

These ideas, and their application to eikonal models of hadronic scattering (built out 
of eikonal models for the underlying quarks), should be of use to phenomenologists and 
experimentalists, who must translate pure QCD theory into practical predictions and require 
a separate analysis of binding, and that topic is not covered in this paper. Rather, we confine 
ourselves to the basic properties of quark scattering by the multiple gluon exchanges noted 
above, and observe qualitative results depending upon the impact parameter, which are 
reminiscent of the MIT Bag Model, and of asymptotic freedom. 

One of the common features of QED is that MGI is incompatible with manifest Lorentz 
covariance (MLC); that is, one must choose, and has always chosen, a gauge-dependent 
formulation as the price of MLC. The reasons are well known, stemming from the effective 
balancing act of constraints vs. true degrees of photonic freedom. Traditionally, it has been 
most convenient to choose a gauge dependence for the covariant photon propagator, assure 
oneself of the gauge independence of radiative corrections to that photon propagator (by 
means of rigorous fermion-charge conservation), and accept the necessity of gauge-dependent 
photon propagators as long as all properly-defined S-matrix elements of the theory can be 
shown to be independent of gauge \L\. In QED, as Schwinger has shown [s]. Green's functions 
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of operators calculated in the Coulomb gauge can be transformed into Green's functions in 
conventional relativistic gauges by adjoining an operator gauge transformation to the original 
operator, and so retain the basic quantum formalism without the need for indefinite metric 
quantization. 

The same, basic formalism may be followed in analytic treatments of perturbative QCD. 
An additional, additive feature has been the apparent necessity of adjoining "ghost" fields 
to the theory, in order to produce a conventional representation of the gluon propagator 
with its proper degrees of freedom 1^. As shown in this note, there exists a simple, 'virtual- 
gluon-ghost' interpretation which can be used to 'spark' a MGI and MLC formulation of 
quark/ ant i- quark interactions; and in this formulation, the non-perturbative, mathematical 
representation of physical processes is expressed by an FI over the position and color co- 
ordinates of a single, anti-symmetric color tensor, x1j.u{^)- This integral, which, long ago, 
was suggested by Halpern |9|, [lO|, begins life in the definition of an FI; but due to the ghost 
nature of the virtual gluons, is reduced to a single n-fold integral over 'local' position and 
color coordinates. One advantage of the present method is that it is accessible to couplings 
of any size; and, in fact, the calculations appear to simplify considerably in the limit of 
strong coupling. 

As possibly the simplest, non-trivial illustration, we set up the calculation of a high-energy 
Q and/or Q eikonal scattering amplitude, in quenched approximation, and at large coupling, 
using recent eikonal techniques for non-Abelian interactions [2]. At the end of this model 
calculation, one can see that Halpern's integral describes an effective, 'almost-contact' in- 
teraction between the Q's and/or Q's, replacing the conventional, boson- propaga tor 'action- 



at- a- distance' Abelian eikonal result. As in QED, or any Abelian theory Il2l. Il3l. \l4. Il5|. 
a logarithmic growth of a total cross section will require at least a partial lifting of the 
quenched approximation; in this simplest model invoking quenching, one finds a scattering 
amplitude dependent only upon momentum transfer (or impact parameter), with reasonable 
color structure. 

One simplification employed below should be stressed, for although the model we present 
resembles an eikonal calculation, certain rather complicated normalization factors have, for 
convenience, been omitted, as noted at the appropriate place. The thrust of this presenta- 
tion is therefore limited to display the method of virtual-gluon exchanges; and to show, in a 
scattering context, how dependence upon momentum transfer or impact parameter controls 
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the disruptive effects of color fluctuations on otherwise-coherent, eikonal-like exchanges. To 
put this calculation into a strict eikonal framework, as in Appendix B of the QA refer- 
ence jsl, one need calculate the neglected normalization factors; and hence when we refer 
to the 'scattering amplitude' we mean an unrenormalized, MGI and MLC quantity whose 
magnitude we can only compare for different impact parameters. 

A list of abbreviations of frequently-used phrases has been added as Appendix [Dl 



II. FORMULATION 



Begin with QED, and its free-photon Lagrangian, 

C, = ~-jl = --id,A,-d,A,f. (2) 



Its action integral may be rewritten as 



^'xCo = ~\l {d.A,f + \ I {d,A,Y (3) 



11 A, {~d')A, + \j {d,A,)\ 



and the difficulty of maintaining both MGI and MLC appears at this stage. What has 
typically been done since the original days of Fermi, who simply neglected the inconvenient 
{d^A^)^ term, is to use the latter to define a relativistic gauge in which all calculations 
retain MLC, while relying upon strict charge conservation to maintain an effective gauge 
invariance of the theory. 

The choice of relativistic gauge can be arranged in various ways; and what shall be done 
here, in the context of the preceding paragraphs, is to multiply this inconvenient term by 
the real parameter A, and transfer it into an effective 'interaction' term. For definiteness, 
begin with the free-field, (A = 0, Feynman) propagator Tifjiu = Sfj^^Dc, where (— 9^)Dc = 1, 
and the free-field Generating Functional (GF) 

2(o)^]=,i/rf..^ (4) 

and operate upon it by the 'interaction' A-term, to produce a new, free-field GF 

zPlj] = e'^'f^'^^^^' .eiS^^'-= (5) 



A- 



.1 A 

i Sj 
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/i-D^^.j -|Trln[l-Af] 



where Dcf^^ = {6^1, — (d^du/d'^) Dc, with ^ = A/(A — 1). The functional operation of ([5]) 
is fully equivalent to a bosonic, gaussian, FI; and such 'linkage operation' statements are 
frequently more convenient than the standard FI representations, since they do not require 
specification of infinite normalization constants. 

The Tr-Log term is an infinite phase factor, representing the sum of the vacuum energies 
generated by longitudinal and time-like photons, with a weight A arbitrarily inserted; this 
quantity could have been removed by an appropriate version of normal ordering, but can 
more simply be absorbed into an overall normalization constant. 

Again starting from the Dc°/ljy of a Feynman propagator, and including the usual fermion 
interaction Cmt = igi'l ■ ^i^-, ^^^d the gauge 'interaction' |A it is also easy to show 

that one generates the standard, Schwinger functional solution in the gauge 



-,(0) 

■37 

where the phase factor of ([5]) has been absorbed into M . It will be convenient to rearrange 
dn]) using the easily-proven identity 



i 6j 



where = -| / jx' Dc^^ ■ ^, so that % now reads 



^h^-^ ^ . e^i^' . T\A\ ^ , (7) 



This is the functional QED we know, and have used for a half-century. 
We now come to QCD, with 

£ = - JfJ, - ^ [m + 7 ■ 9 - ^^77 ■ ^ ■ A] V-, (9) 

and F;J^ = d^Al - d^A"^ + gf^^^A^Al = f;j^ + gf'''"'A^^Al. Since 'proper' quantization in 
the Coulomb gauge, for the free and interacting theories yield the same ETCRs for QCD as 
for QED (with an extra 6ab color factor appearing in all relevant equations); and since at 
g = 0, QCD is the same free- field theory as QED (except for additional color indices); and 
since QED in any of the conventional relativistic gauges can be obtained by treating the 
I A J (df^Af^)'^ as an 'interaction' (as above); and therefore, rather than re-invent the wheel, 
we set up QCD in the form used above for QED. 
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As a final preliminary step, we write 




with f;, = d^Al - d,A1 and 

C'[A] = — ^(2f^"^ + gf"-'^'^A^^Al){gf"-^^A^Al); and for subsequent usage, after an integration- 
by-parts, we note the exact relation 

- i = /^^ (-^^) + {d.A'^y + 1 c'[A], (11) 

[In the next few paragraphs, for simplicity, we suppress the quark variables, which will be 
re-inserted at the end of this gluon argument.] 

To choose a particular relativistic gauge, multiply the 2nd RHS term of (ITTl) by A, and 
include this term as part of the interaction, to obtain the familiar QCD generating functional 
(GF) in the relativistic gauge specified hj ( = X/{X — 1) 

zi^Uj] = ATe^^^'t^^] ■ e^'^^'^'^^ . ei/^-°^ (12) 
or after rearrangement 

2SDL7]=Are^^-'[^*]-e^/-°^^'^ (13) 

with the determinantal phase factor of ([5]) included in the normalization Af, and a 6°''^ 
associated with each free-gluon propagator oif^l^*. 

After re-inserting the quark variables, and after rearrangement, expansion of (fT3|l in 
powers of g generates the conventional Feynman graphs of perturbation theory in the gauge 
(. If one wishes to have a conventional form for the renormalized gluon propagators, one 
can insert Faddeev-Popov ghosts into the Lagrangian. But it is clear that all choices of 
A are possible except A = 1, for that choice leads to C ^ cxd and an undefined gluon 
propagator. This is an unfortunate situation because the choice A = 1 is precisely the one 
which corresponds to MGI in QCD, as is clear from ffTTl) . 

But there is a very simple, alternate way of writing (fT3|) . by replacing the J C of that 
equation by the relation given by the exact ( fTTll . 
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which (continuing to suppress the quark variables) yields 



(15) 



i Sj 



It is now obvious that the choice A = 1 can be made. It will become clear below that, in 
the form f|T5|) . these operations are exactly equivalent to the introduction of a gluonic ghost 
field; and it is this 'ghost property' for virtual gluon exchanges that generates an exceedingly 
simple, MGI and MLC result for the present eikonal model — and for all subsequent radiative 
corrections to this model that can be written. This ghost mechanism occurs because the 
ghost gluon has been introduced by the Feynman propagator assumption which leads to 



the factor exp 



(0) 



J 



Mj lA^iA is that 



of ([IS]), while the term exp [| / A"^ (-8^) A"^] 
functional operator which will remove every such propagator from the sum of all virtual 
processes of every n-point function of the theory, without exception. In effect, the gluon 
ghost acts as a 'spark plug' to generate the MGI and MLC interactions of the theory, which 
then take on a remarkably simple form. 

If one argues that because no color gluons can ever be asymptotic, it is then reasonable to 



»(0) 



of the rearranged OF; or, if one wishes 



suppress the leading RHS factor exp | / j • Dc ■ J 
to retain the specification of individual gluons, that factor may be retained, and standard 
Faddeev-Popov ghosts inserted to guarantee its proper perturbative renormalization. In the 
example to be given shortly, this factor plays no role and will therefore be omitted. 
After rearrangement, and after re- inserting the quark variables, f|T5|) becomes 



A=/D 



(o).'(16) 



and we next invoke the representation suggested by 
Halpern 9|, 



e 4 



f^'=Af' y dMe^^^^^-^'+s/^M.FS., (17) 

where / d[x] = lit Ila n^>jy / ^X^i'^Ji), so that (IT7|) represents a functional integral over 
the anti-symmetric tensor X^jy('W^)- Here, all space-time is broken up into small regions of 
size 5^ about each point wi and M' is a normalization constant so chosen that the RHS 
of (fT7|) becomes equal to unity as FJ^^ — > 0. In this way, the GF may be rewritten as 
(AT' • AT = U" M) 



i J f)-Gc[A]-q+\.[A] 



(18) 



■ e 



where S)r = -t/^.Di°).^. 

As noted above, we treat the quarks and anti-quarks as stable entities during the scatter- 
ing; and then must calculate functional derivatives with respect to the sources f]{xi), i]{yi), 
fj{x2), and ^7(2/2), which bring down factors of G\{xi,yi\A) and Gj(x2, ?/2|^), where the su- 
perscripts I and E refer to the scattering fermions. With standard mass-shell amputation, 
we pass to the small-momentum-transfer limit of the eikonal model [sl, derived in detail in 
Appendix B of this reference for the specific case of QQ scattering, and using the conven- 
tional, FI approach in an axial gauge. (The discussion of Appendix B of ref. jsl contains the 
full QCD, with cubic and quartic gluon interactions.) 

The quark scattering amplitude is given by the familiar eikonal form [3|], 

T(.,t) = ^ |d^6e-^[l-e«(-^)], (19) 

s = -(pi +P2)^ 

. / i\2 2 CM ^2 

t = -[p^-p^) = -q > -q , 

while the exponential of the eikonal function, E = exp [zX], is obtained in this quenched for- 
malism by the appropriately normalized (as in (B.32) of this reference) action of the linkage 
operator: exp [-f / ^ ■ 0^°^ ■ ^] upon exp [f / x ■ F + f / (-0^) A^l] ■OE{p,,p[,p2,p'2} 
in the limit A ^ 0, where the factors denoted by 0E{- ■ ■ } are the ordered exponentials con- 
tributed by the Green's functions corresponding to the incident and outgoing particles, as 
noted below. As remarked in the previous section, for simplicity of presentation, certain 
normalization factors shall be suppressed; and therefore our result is only a qualitative ex- 
pression of the scattering amplitude, or rather, of the eikonal exponential E = exp[«X], in its 
dependence upon impact parameter. But from this qualitative result it will be possible — by 
means of two numerical integrations — to obtain a qualitative picture of the effective inter- 
action potential between a pair of quarks or of a quark and an anti-quark. 

In QED each such Green's function Gc[A] would contribute an exponential factor 
exp [ig f d^w7lf^{w)A^{w)~\ with 

/O poo 
ds6{w - y + sp) + p'^ / ds6{w-y + sp') (20) 
-00 Jo 

/+00 
ds 6{w — y + sp), 
-00 
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but in QCD it will generate an ordered exponential (OE) of form 

(^exp igp^ dsA';^iy-sp)X'^^ . (21) 

In order to extract the A-dependence from such OE, we rewrite (1211) as 

d[a]S[a'^{s) - gp.A^iy - sp)] ■ (^e^/.+^dsA'^.n.))^ , (22) 

where / d[a] is a functional integral defined over all values of the mesh coordinates — oo < 
Si < +00. One then writes a representation for the (5-functional of (|22|) . so that the OE of 
( |2T|) becomes 



■'Jd[a] Jd[n] 



-i dsn-is)[a-{s)-gp^A-{y-sp)] _ f i ds \''a-(s) 



(23) 



where J\f' = ^ normalization constant for the functional integral over the 'proper 

time' s values, with the width of each mesh given by A (of dimension (length)^), A ~ 5^. 
These operations have become routine in eikonal analysis 

In our present QCD eikonal scattering amplitude, each Q, or Q, is described by a Green's 
function Gc(xj, yi\A), and each has an OE of the form expressed by ( 123|) . with corresponding 
Pi, yi, and ^^^(sj) variables. In QCD eikonal models, only the interaction corresponding 
to multiple gluon exchanges between the scattering QQ are retained, and the functions 
contributing to the eikonal amplitude will contain those pair-wise-interaction variables in 
the manner of (!36|) . below. 

There is another A-dependence contribution to Gc[A], the OE denoted by 



(exp g J ds a^^ F^^^iy - sp) ^ 



(24) 



where this OE is again defined by its s-value. However, in the present virtual-ghost-gluon 
calculation, a simple scaling argument shows that these spin-sensitive terms do not appear. 

We retain the basic idea of an eikonal model, concerned with the interaction of a Q, 
or Q, each treated as a particle of renormalized mass m and color charge g; and to effect 
this statement, we suppress all self-energy structure of each Q, or Q. We also suppress the 
L[y4]-dependence, as in a quenched approximation, where the scattering is assumed to occur 
so quickly that charge renormalization effects and any change in the fundamental vacuum 
structure have insufficient time to react. As shown in the original calculation of Cheng and 



Wu 



11 



12 



13 



15|, contributions from Li[A] are essential for the increase of total cross 
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sections with scattering energy; and such effects will be missing in the simple model here 
described. 

From (1181) and (1231) . our eikonal exponential function E = e*-^ will be proportional to 



d[x] 



3 4 / X 



S r\{0) s 



-1 f -2-.D 

e 2 J «A SA ■ Q2 



where 



if ■ x)t = r'^x' 



The linkage operation is again Gaussian, and yields 



-iTrln Tl-^-D™^ 



■ e 



D 



(0) 



-,(0) 



■Q 



and here is where the 'ghost-magic' appears, since 



1-/C-D(°) = l-[^7/-X+M')]Di°) 
= -^(/■x)DW. 



The Q-dependence of ( l26i) is then just 



2^? 



Q-{f-xr-Q 



with all Yif^ propagators canceling away, leaving an integral over a single space-time variable, 

-Ygj^""" Q» [/ ■ ^(^)rxt (29) 

Precisely this form of effective interaction was previously found in an instanton approxima- 
tion to a QCD field-strength formalism 6| with the difference that our fl29|) does not contain 
a gauge-fixing term of form ^[(^"(a;)] of Eq. (17) of that paper (using the notation of that 
paper, where the original gauge-fixing dependence of the variables was replaced by a 



(25) 



(26) 



(27) 



(28) 
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simpler gauge fixing of tlie F^^,, and tlien transferred to the field). We consider this an 
unphysical difference because of the following argument. 

Different gauges are traditionally chosen in order to simplify calculations of different 
processes; for example, in the calculation of infrared effects in QED, the Yennie gauge, 
( = —2, is to be preferred over that of any other ( value because it simplifies the analysis. 
Gauge A is chosen to calculate process A because gauge B or C would entail a great deal of 
unnecessary work; but because the Physics is independent of gauge, if no errors are made, 
then the use of gauge B or C or any other gauge must give exactly the same results as does 
the use of gauge A. Instead of using gauge A to calculate any process, we could use the 
average of gauge A and of gauge B, or take the average over all possible gauges; the physical 
answer must be the same. In the language of reference 6|, the choice of g"'{x) is arbitrary; 
and so we suppose that we may average over an arbitrary number of such g°'{x)] and if there 
are a continuous number of such functions, as can surely be imagined and constructed, then 
instead of using a particular g°'{x), we may simply calculate / d[g] 6[g{x)] over the complete 
functional space of such functions, and divide by the (infinite) volume of such a space, which 
latter quantity may be absorbed into an overall normalization constant. The result of this 
last summation, taken under the / d[x] integrals, is just unity, which is the form of our fl29l) . 

Just as an average over all paths is path-independent, so an average over all possible gauge 
choices is gauge independent; our result is gauge 

independence was forced by the ghost mechanism, which automatically removes the gluon 
propagators carrying any arbitrary choice of initial gauge. The restriction initially made for 
the use of the Feynman gauge for Dc°'* was only for simplicity of presentation; the entire 
discussion could have been carried through by adding and subtracting another gauge-fixing 
term to the Lagrangian, adding it to Cq and subtracting it from C. Again expressing the 
interaction in terms of Halpern's integral, one finds exactly the same cancelations, except 
that it is a propagator in an arbitrary gauge that is removed; the details are in Appendix O 
At this point it may be useful to digress into just how ghost enter QFT, especially in the 
context of a linkage operation. For immediate relevance, consider a bosonic ghost, which 
might be introduced in order to have an intuitive representation of a determinantal factor 
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exp [— |Trlni3], where B is any desired quantity, or operator. Consider the operation 



iTrln(l-Dt,X) 



(30) 



with the choice /C = ^ + i3. Then, because Dc ■ ^ = 1, 1 — Dc ■ /C = — ■ B, and 
Dc(--B ■ Bc)-^ = -B-\ so that ([30]) becomes 



mahzation constant, the remainder is just the desired term. Of course, fl3Tl) is completely 
equivalent to the Gaussian functional integral 



with an appropriate normalization A/". But the ghost mechanism is not simply just the 
Gaussian integral fl32l) . for in the linkage operator formalism one sees the removal of all Dc 
from the internal, virtual structure of the theory. It is in this sense that our virtual QCD 
formalism corresponds to the use of a gluon ghost, which has been forced upon us by the 
requirement of MGI; and that requirement is then rigorously satisfied by the removal of all 
gauge-dependent gluon propagators. 

It should be emphasized that this ghost removal will occur automatically for every cor- 
rection, quenched or unquenched, to the simplified limits of this example. For example, 
the i^lA] terms neglected in this quenched calculation can be retained by a straightforward 
expansion of exp{L[y4]} in powers of L's; and every L[A] so included may be expressed in 
terms of an exact Fradkin representation [sl, , which is itself not more complicated than 
a sequence of operations upon an exponential of linear and quadratic ^-dependence. The 
totality of such radiative corrections, exactly or in any form of approximation, will always 
retain the same form as in fl29|) or fl3T]) above, with JC and Q having added terms; but the 
removal of all Dc propagators must again occur, as MGI is maintained. 

The question then arises: If the gluon propagators are to disappear, what is going to 
replace them as the 'carriers' of interactions from one Q, or Q, to another? And the an- 
swer is that the Halpern field x'Jj.u takes on a new and physical significance as the carrier of 
the totality of virtual-gluon interactions, in the form (as seen below) of an 'almost contact' 



Setting A = 0, and treating exp [— |Tr In (— Dc)] as an unimportant, if divergent, nor- 




(32) 
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interaction. And it is this novel interpretation which has the abihty, in one succinct if comph- 
cated representation, to display the effective QCD interaction for all values of the coupling, 
large or small. To the associated question of a possible phase change for large values of the 
coupling, we, in this paper, make no prediction, for the answer to that question demands an 
evaluation of our results for small g, and comparison with simple QCD perturbation theory. 
Our final Halpern integral is much simpler to evaluate for large coupling, rather than small; 
and in the interests of simplicity of presentation, that question has been left unanswered. 
Another untouched question is whether color-charge renormalization in this formalism will 
require additional Faddeev-Popov ghosts, which could certainly be inserted, if desired. 
Returning to (123]) and fl26p . now written in the form 

U j d[x] e^^^ ■ [dei{gf ■ x)T'^ ■ e'i ^ S''^ '^]"'-®. (33) 

We first drop the self-energy parts of the Q-dependence, retaining for the exponential factor 

-^9' I'J^l-igf-xr'C-'^l (34) 

and then, for simplicity, discard all but the largest ^f-dependence of (IMl) . 

-^9' jni-{gf-x)-X.-^'^^- (35) 

Inserting the eikonal representations of IZ^^ and T^j,^, in the CM of the scattering quarks, 
we need to evaluate 

dsi / ds2Pi,P2u ■ n^ds,) ■ / d^^ [f-x{w)r'\l, ■ n'^{s2) (36) 

OO J —OO J 

xS^^\w -yi + sipi) • 5^'^\w -y2 + 52^2)- 

Here, pi, p2, yi and y2 are the relevant 4- momenta and space-time coordinates appearing in 
each Q/Q Green's function, and they are evaluated in the CM of the scattering Q/Q, which 
are (initially) assumed to have zero relative transverse momentum. In this way, 

Pl,4 = P2,4 = iE, pi,3 = -P2,3 = P, (37) 



Pi,i = Pi,2 = P2,i = P2,2 = 0, z = yi-y2 

2 p ^ E 



if -X) ^|34 = i{f -X) ^IsQ, s± = - 
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so that the product of the two delta-functions of (136|) becomes 



5^^\yi,L - y2,±) ■ S{si - s+) ■ 6{s2 - s_) (38) 
■6^'^\w± -y±) -6 (^wl - ^ {yi,L + 1/2,l) j 



E \ 1 

■0 [Wq- yin + — 2/1, L 



where y_[_ = yi^± = —y2,i. = and the zero of CM time is chosen when both particles 
are at their distance of closest approach, when i/i q = t/2,o = 0; then, for all times, zq = 
2/1,0 — 2/2,0 = 0. Hence, si = S2] and since 2/1,0 = •ymsi, si = 2/1,0/(7^^)5 and for large 7 
and any reasonable duration of the scattering, si ~ ~ S2. Also, 2/1,0 + 2/2,0 = 22/o, and 
2/1, L + 2/2, L = 0, and the entire fl5^ may be written as 

tg6^'\b)nm [f-xiw)rXo^UO), (39) 

where the expected anti-symmetry of the /i, u variables of [/ ■ x]^^ has been used, together 
with the pi,^, p2,u values appropriate to the CM. Note that the w variable of [/ ■ x]~^ is a 
fixed 4-vector, given by w^^-* = {y±, 0^; 2/0) for E/p ^ 1. 

This last restriction immediately means that only this w^^\ of all the possible w- values of 
the original Halpern representation, is relevant to this interaction; and all of the other w^- 
terms of that functional integral, with their normalization factors, are effectively removed 
from the computation in the form of an uninteresting, convergent, normalization factor, 

Af' I d[x] ei^^ sjdei{gf-xr\ (40) 

which separates itself from the 6-dependent part of the calculation. The latter, in contrast, 
is now given by 

n n [tA^I r Vdet(,/.x)-^et^^ (41) 

^ gi3<5(2)(fe) Uf (0) If-xMr' I ll f^l (0) ^ 

where A = 5^, and 6 refers to the small distance, in each of four space-time directions 
surrounding the point w^^\ For a later purpose, we shall borrow a divergent factor from 
one of the normalization terms, and so insert a A, multiplying x, into the determinant of 
(|4T1) . Henceforth, we suppress the w^^^ symbol, with the understanding that the integral of 
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fHT]) refers to the summation over all possible values of the quantity we write 

the measure of (jUJ) as Y[ij.>i/ 1 ^"'Xfj.u, where n refers to the number of independent color 
contributions of SU(N). 

The next step is to rescale the A-dependence of fj^Tl) . defining xfiu — '^X% ^^'^ ^o obtain 

(2vr)-2 n / d"^/^- 7det(^/ ■ (42) 

and we must then interpret the quantity For this, first write a Fourier representa- 

tion 

5(2) (6) = i2Tifj d^k±e'^^\ (43) 

and realize that this integral requires a specification of all k±. But is this reasonable in an 
eikonal model of quarks, where we understand that such quarks can never be measured in 
isolation, with precise values of momenta? Rather, we must extend this eikonal model to 
allow for unmeasurable transverse momenta exchanged between quarks of the same hadron, 
before any quarks in different hadrons can be imagined to interact with each other, as 
is the conceptual situation of this calculation. That transverse momenta, which can be 
treated as an average quantity even though it can never be measured with precision, will 
certainly be smaller than the CM momenta, or the CM energy, of the hadrons which are 
actually scattering; and it will be on the same order of magnitude as the transverse momenta 
defining the (5-function above. In other words, taking into account that we are talking about 
quark scattering, rather than particle scattering, the magnitude of the transverse momenta 
inside J d'^k± must be limited; and the natural parameter which sets the scale for high- 
energy scattering, in which eikonal models are most relevant, is the CM scattering energy 
of the hadrons. We therefore insert under the integral of / d'^k± a limiting factor for its 
transverse momenta; this can be done in many, physically-equivalent ways, but perhaps the 



simplest is to use exp 



-kl/M^^ 



with M on the order of the CM scattering energy. This 

, and has 



replaces 5^'^\b) by a more realistic Gaussian distribution (M'^/Atc) ■ exp — M^6^/4 
the further advantage that the product AS^'^\b) is now proportional to the dimensionless 
quantity AM^. Since it was our eikonal model that, in part, defined A, it is reasonable to 
choose the product AM^ = ^ as a number ~ 0(1), thereby replacing the original A5^'^\b) 
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by Lp{h) = ^ exp |^-M262/4j . Eq. (02]) then becomes 

(27r)-2 n / \/det{gf-xr' ■ ei^'+^^^^^^^^^o) [f-xn:>i(^). (44) 

III. EVALUATION 

We next turn to the evaluation of fj44|) . which is to be inserted under the normahzed 
functional integrals 



Af" j d[ni] J d[nji] exp 



—I 



— oo 



(45) 



But the 6-dependence of f H4|) is associated with fii(O) and i^i(O); and this means that all of 
the other Si values, Sj 7^ 0, of the / d[Qi] and / d[Qj\ may be integrated immediately, with all 
yielding factors of 6{ai{si)) and 5(an(sj)), 7^ 7^ Sj. Only the normalized contributions 
of each functional integral with s = are relevant here. In the Quasi-Abelian model of 



reference js], this was suggested as an intelligent approximation for the SU(2) eikonal model 
considered there; but here, for SU(3), it is an almost automatic consequence of the 'locality' 
of the gluon-ghost mechanism. And it has the extremely convenient effect of transforming 
the remaining OE integrations over 

(exp [—i J «! • Ai] ) ^ and (exp [—i J aj ■ An] ) ^ into ordinary integrals over unordered quan- 
tities, 

j d"«ie-'^°i(°)-^i (46) 



271^ 

with the result that all that remains of the color dynamics are the tedious but straightfor- 
ward, ordinary integrals 

(27r)-2" J (Tai J d"an j d'Tli j d"l]i (47) 

where we have rescaled the ai,n variables, and suppressed their now useless (0) notation, as 
well as the notational change: X ~^ for the remaining x integration. 
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To evaluate (H7|) one needs a representation of the inverse of the doubly anti-symmetric 

I ab 

3n(n — 1) independent quantities comprising this quantity, and the simplest, compatible 



matrix, [/ ■ x] ^ | • If there exist n color and 4 space-time coordinates, there are then 



assumption is to write 



where we expect G"''' and H^u each to be anti-symmetric. If fHHj) is true, there then follows 
the necessary condition 

<5%. = $^(/-x)irAG'^'-^A., (49) 

c,A 

which can be used to provide implicit representations for both G and H, as follows. Set 
a = 6 in fHOD and sum over all color coordinates to obtain 



■ Hxu, (50) 



X 

from which it follows that 



6,c 



Similarly, set /i = z/ and sum over space-time indices to obtain 

{G-T = \j:(f-X)CH>^.- (52) 

fj,,\ 

It will be convenient to define 

= - E f'^'G^'xlu ^ E ^'^M- (53) 

b.c e 

SO that and its inverse can be expressed as 

{H-'),^ = Q,^. H,. = {Q-')^^. (54) 

The general statement of the inverse of an anti-symmetric, 4x4 matrix can be used to 
represent Q^^ as 

but it will be most useful to note that only one of the six, independent H^j^^, H^o, is multiplied 
by the factor gip in the exponential of (jHj); and for small b and large g, this contribution 
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will be large. Does this carry the implication that, for all color indices, the Xso will typically 
be larger than the Xa/s '^^ other Lorentz indices? Not necessarily, but in the interest of 
simplifying the computations we shall assume that only X30 and Xu are of interest. [In order 
to prevent det[(5] from vanishing, it is necessary to retain one other Xap in addition to Xso-] 
With this approximation, det[(5] — * QuQlo^ 

= QMi ■ 



SO that 



H3o={Q-'),, = ^, (57) 

3 



^30 

and 

with G"''^ given by the inverse of fl58|) . 

Does the inverse of exist? The inverse of an anti-symmetric matrix M"* of eight 
rows and columns is given by 

[M-^y'' = ^[detM]-5 ^ahcdefgh^cd^ef^gh^ ^gg^ 

where ("■^'^'^^fah jg wmi anti-symmetric tensor of eight dimensions. However, if M°'^ = 
jabc yc ^ ^j- . yyb^ calculation shows Q that det[/ ■V] = 0, for any and every value of 
the color vector V^. In general, inverses of such Lie- valued sums do not exist, and it might 
appear that this MGI calculation must grind to a halt. However, what is relevant is the 
combination of G with H, not G alone; and therefore let us give a physicist's redefinition of 
the problem. We shall define the determinant of a matrix M as: det[M] + A^, where A — >■ 
as a subsequent condition. It is also understood that the elements of M are dimensionless. 

Rewriting ( l52l) in the form (G~^)"^ = (/ ■ V)""^, where = jYj,_ixXI.x Hx^i, the corre- 
sponding G°'^ may be expressed as G"^ = G°'^/X, where G"''' is defined by (1/48) multiplying 



the corresponding numerator of (!59|l . with M = (/ ■ V) . From (l53l) the quantities Q 
and may be written as Q^j^u = Q^iu/X, (f = (t l^-, where Q and q are defined in terms of 
the finite G. Then, from fl54p one may write if^j, = A (0^^)^^ = ^^H^^] and in this way the 
product G""^ -H^u of fHHl) becomes G"''^- H^^, and is independent of A. Only G"^ quantities are 
needed in the subsequent analysis of color dynamics; although one finds a factor of det[G~^] 
required at one point in the calculation, it is immediately followed by a factor det[G]; and 
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the product of two such determinants is unity. But were there a divergent contribution of 
any form associated with the original (/ ■ which we have represented by the product 

from G ■ H, there exist separate arguments to show that such divergences have no^ effect 
on the Physics; one of those arguments appears in the paper by Reinhardt, et al. [6], and 
an independent proof is given in Appendix [X] of the present paper. Perhaps the simplest 
argument is to observe that any singularity of (/ ■ x)~^ will cause the exponent of (14^ to 
oscillate infinitely rapidly, and make no contribution to the integral. 

Finally, in the special limit of small impact parameter, the eikonal exponential E[x] will 
reduce to a finite set of possible terms — all considerably smaller than that of the large impact 
parameter result — involving the magnitude of the diagonalized components of G. In this 
way, because the color coordinates are coupled to space-time, the procedure is well defined 
and yields qualitative results in agreement with QCD intuition. We shall find that for 
large impact parameters the scattering is coherent, with the quarks retaining their original 
color, while for smaller impact parameters, color fluctuations reduce the magnitude of the 
amplitude. 

IV. ESTIMATION 

If the Halpern variable x% is written as 3^^^ r^^, where 3^^, represents the color-projection 
of a 'magnitude' r^^, inspection of the original inter-relations of G and H shows that G 
is independent of the 'magnitudes' r, and depends only upon the 3; and we shall assume 
the same dependence for H and G. In contrast, while dependent upon the 3, varies as 
the inverse of the r variables; and it is this latter r-dependence which appears to be most 
relevant to the overall color properties of the amplitude. We shall therefore treat the 3- 
dependence as producing relatively unimportant averages which are to be relegated to later 
numerical integrations, and concentrate in what follows on the output of the r-integrals. 
Since the ^^(yj-dependence of f HTl) is associated with the dependence of i^so, integration over 
X12 variables can be moved into a separate, uninteresting normalization constant; and we 
suppress the (30)-subscripts of the remaining X30 variables. 

The (7y9-dependent exponential factor of (l47l) is then 

e^^[ig^Q.tG''^Q!'Jr\ (60) 
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and we first consider the integral 

{2rr)-" J d«fiie-*"i^?e^^^^?'5"^^/'- (61) 

Now define fi^ = {0-^)^"%, so that dHID becomes 



and / d"'fli yields 



('i_Vdet[G-']e-<'i(-'-')"'°S(lfe) (63) 

Since ^ i [(^12)^ ~ (Xso)^]' after removing the Xi2"dependence, there remain the 
(7V9-dependent integrals 

(2.)-»/d.V,e-".-../a"„.e-.Vj. (64) 

xdetfG-i] ■ j d"xe-*"'/' (^^^ 



where X30 = = rf, 



/r+oo 
d"x = n / (65) 

/ + 00 POO ^ 

,. -OD Jo 



and with dx^ = rd^'^, 

"+1 



(67) 



d"x^2ny^ d3'=5(l-^(3'^)2) drr"-\ (66) 
Then, ( 16^ may be rewritten as 

c -^^1 a 

xdetfG-^] ^ drr"-i (^^^ e"'"'/^ 

x/d"«.e— 7d"«.e— ^ 
xe-<^)"S('5-)"^?. 
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For clarity of presentation, in the passage from (jUj) and (HSil to (162|) . we have suppressed 
the factor of a/ det{gf ■ x)~^ of fj44l) . From the discussion of Section [TTll and that of Appendix 
Rl this omitted term will contribute a factor of r^^/^ to the integrand of (1671) . which will 
have no bearing on the qualitative conclusions of Sections [IV] and |Vl 

It will now be most convenient to isolate the ai^n factors from the Ai.n factors, by writing 

g-iAi-ai _ (271^ j (Tv j d"fie*^ (''-"i) •e-*''-^! (68) 



so that integration over the a\j may be performed, 

j d"ai y d"an e-^^^'^"^"^-^ ■ g-^^K'^-O"'"!/" (69) 
= j d"aie-*""-^5(afi° - (G-^)"^4), 

where a = gf/r. With the variable change: = G'"^j3'^, this becomes 

(27r)Vdet[(5] y"d"/3e-^^-'5-^5(/?-afi) = (27r)Ve-'"^-'5-^, (70) 



and one notes that the determinantal factor of f[70|) combines with that of fl67j) to produce 
a factor of unity. 

At this point is useful to perform the remaining v, w integrals written in the form 

J (Tn j (Tne-'''^-^-^ Ji{n)Mn), (71) 



where 



Ji(fi) = (27r)-" j (Tve-'^-^' ■ e'''-^, (72) 
and 

Ji(n) = (27r)-" j (Tw e-'""-^^ ■ e™'^. (73) 

Clearly, for g^{h) 0, (!7T!) reduces to a constant, independent of color factors, so that in 
this limit the initial and final quark colors must remain the same; but for large gf{b), there 
will be oscillations involving changing color coordinates away from that constant, so that 
the magnitude of the 6-dependent amplitude will be reduced. 

We have carried out a simple estimation of this effect for the simplest case of SU(2) in 
Appendix [Bl and find that the expectations described in the above paragraph hold true: 
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Color fluctuations at small impact parameter diminish the coherent scattering produced at 
larger impact parameter. This non-perturbative and gauge-invariant statement can form the 
conceptual basis of quark scattering and binding, as well as asymptotic freedom. A more 
precise statement must await a careful program of numerical integration, which we are not 
able to perform. But there can be little doubt of the qualitative nature of the output of 
such a detailed calculation; and for this reason, we believe that the methods described in 
this paper open a door to the realistic estimation and calculation of detailed QCD processes, 
properly gauge invariant, and containing all orders of coupling. 

V. SUMMARY AND EXPECTATIONS 

The above Sections have described a new method of calculating a particular scattering 
process in QCD, to all orders of the coupling and with CI and LC assured. We have 
made a number of approximations for ease of presentation, as well as for our inability of 
performing certain relatively unimportant integrations which must be left for subsequent 
numerical integration. Our result is a qualitative expression of the eikonal exponential 
function E = exp [iX] , given as a function of the square of the impact parameter between the 
scattering particles. And from this quantity, by a process requiring numerical integrations, 
it is, in principle, possible to obtain a qualitative idea of the effective interaction potential 
between a pair of quarks or of a quark and an anti-quark. 

To see this, return for a moment to the simple, potential theory problem of a particle 
scattering from an external potential T^(|r|). There, the corresponding function E is given by 
the exponential of a simple kinematical factor multiplying the two-dimensional expression 
of that potential, obtained — as a result of a relevant, eikonal-calculation prescription — by 
calculating the three-dimensional Fourier transform of that potential, and setting 

the longitudinal component of that 3-momentum equal to zero, to obtain V^(|/c_l|). 

In all previous field theory models, or approximate calculations of subsets of Feynman 
graphs, which yield eikonals, x(&), dependent upon the square of the impact parameter, the 
two-dimensional Fourier transform of that eikonal generates an effective V^(|A;x|); and the 
simple 'extension' of |A;_l| to the full, three-dimensional \k\, produces the Fourier transform 
of the original potential y(|r|). The same process may be considered for the log of the 
function E we have obtained, with its built-in, qualitative approximations. Because of 



24 



the relative complexity of our result, the Fourier transform over its 6-dependence must 
be done numerically; but that is certainly possible, in principle; and it will generate a 
qualitative V^(|A;_l|). Then, the simple enlargement of that argument, from \k±\ to the full 
I A; I, produces y(|A;|); and a subsequent Fourier transform, again performed numerically, will 
yield a qualitative form for the effective potential between quarks and/or anti-quarks. 

Of course, the potential will, in SU(3), involve Gell-Mann color matrices, as in SU(2) it 
involves Pauli matrices; but these can be included, in principle, in a perhaps tedious but 
straightforward way (as in Eq. (3.8) and the following paragraph of reference js]). Improve- 
ments to our qualitative E can surely be made, by numerical integration over the r- and 
3-factors, at different stages. But here is a method of analytically producing a qualitative, 
effective ^(|r|) — as well as an associated scattering amplitude — which includes contributions 
from every single QCD Feynman graph relevant to the process. 

Of course, we have left out, again for simplicity of presentation, those parts of the Physics 
dealing with charge renormalization, and with the production of particles in the scattering 
process, inelastic effects which have such a unitarity importance to a scattering amplitude. 
But, as explained in the text, these effects can be systematically included in the MGI/MLC 
calculations. They may not be able to be calculated exactly, but it will surely be possible 
to understand their qualitative features. 

The qualitative results seen above for the scattering amplitude — coherent, multiple gluon 
exchange at larger impact parameters, with color fluctuations destroying that coherence at 
smaller distances — are intuitively in agreement with the MIT Bag Model, where quarks are 
'free' when close together but are subject to a confining potential, and tend to bind as they 
move apart; for example, a pion as a bound state of a Q and Q, with the distance between 
them continuously oscillating as they remain bound. Another expected example would be 
the simple vertex function, where the impact parameter of the present calculation becomes 
the conjugate Fourier variable of momentum transfer, so that larger momentum transfers 
correspond to induced color fluctuations and a decrease of the effective coupling strength; 
this is just what would be expected of a theory with 'true' asymptotic freedom, arising from 
the exchanges of an inflnite number of gluons. 

Finally, one must comment on the obvious fact that scattering experiments are performed 
with hadrons, and not with individual quarks; each hadron will involve integrals over the 
transverse momentum or spatial distributions of individual quark wave functions. What we 
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have estimated is the ideahzed case of two quark/anti-quark scattering, suppressing the fact 
that each is bound within its own hadron, and it must be possible to take into account that 
binding. The proper way is to carry out those integrations over the quark coordinates; but 
a simple, physical argument can serve to modify our idealized calculation, as follows. 

Binding suggests that the scale of individual transverse distances, or of the difference 
between those distances of two interacting quarks is controlled by the wave functions, such 
that (b) is never appreciably less than If/i, where /i may be taken as on the order of the 
hadron mass. But our b — B — b, where B refers to the difference of transverse positions of 

— * 

the two hadrons, while b denotes the difference of transverse positions of each quark within 
its hadron. Physically, one expects \B\ > |b|, and the smallest 6- values would be controlled 
by the largest values of j d^k±, which are surely limited, in any eikonal model, by 
the requirement that all transverse momenta associated with, or arising from the exchange 
of gluons must be less than the corresponding longitudinal momenta of the quarks, i.e., 
\k±\<Mr^O(E). 

But there is another question, related to large transverse separations, when the hadron 
amplitude is expected to vanish, because we are fundamentally dealing with short range nu- 
clear forces. How large can the B values become, or how small can the hadronic momentum 
transfer become, before some form of screening sets in and reduces the hadronic amplitude 
to zero? In this case, the needed screening must arise from an interplay of the integrals over 
quark wave functions such that for sufficiently large b, there is effectively no scattering, and 
a 'short-range' force has been achieved; the quark wave functions modify that form of the 
overall, hadronic, eikonal amplitude, such that screening sets in for distances larger than 
1 / n — giving a Yukawa effect between hadrons — while there remains an overall, non-zero and 
coherent scattering for distances less than But from a quark point of view, the essential 
and interesting aspect of our result is that when b becomes so small that b < 1/M, color 
fluctuations begin, and destroy that coherence. 

Finally, one may contrast the qualitative output of such MGI/MLC estimations with 
other, traditional methods of 'summing' Feynman graphs, such as the use of a Bethe- 
Salpcter equation, whose kernel is only known in a low-order perturbative approximation; 
or a renormalization group argument, set up to represent the sum of all perturbative effects, 
but whose beta function is then estimated by a few orders of perturbation theory; or by 
the sum of 'leading-order' perturbative terms, which then omit whole classes of Feynman 
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graphs. In contrast, we believe that the present method holds great hope for generating at 
least qualitative descriptions of field-theory Physics which include, or can be systematically 
made to include, every virtual exchange. 



APPENDIX A: ZERO EIGENVALUES OF [/ x]^^ 

In this appendix, one wishes to get some insight into the role of the operator [/-xl^^'s 
possible zero eigenvalues. One then focuses on the expression (jH]), here rewritten as 

n / dx^odet[(7/-x]~^e^^^°+*^^('')^?[^-^]"l3o^^ (Al) 

a=l ^ 

which is a part of the larger expression 

(27r)-2"y" d"aie-^"i^" j (Taje-'^^'^ (A2) 
X j (Till j d"f]i e-'"'-^' e"*""-^^ 
x(27r)-2j] / dx^odet[^7/-x]-^ 

a 

and where n, as in the main text, is a shortcut for A^^ — 1. One has the relation 

xlo= E (X3o)' = ^tr(x?oA")^ (A3) 

a=l 

\a\h\ _ ATxab 



tr (A^A") = Nd" 

where the A"'s are the n traceless generators of the SU(N) Lie algebra, taken in its n x n- 
dimensional adjoint representation with (A")hc = —if"'^'^. 

Being symmetric under the combined exchange a 6, 3 ^ the operator [/-Xao] can be 
diagonalized and has real eigenvalues. Note that this property applies to [f -Xnv] and can 
be deduced from ( 1331) with Q^, the current given after (p5l) . In the form (lAip . though, this 
property is not transparent. With the p^^^, at i = 1, 2, given after fl36|) . this is because flAip 
results from a re- arrangement of an original expression 

g^{h) X ■ ■ ■ X (pi,3 P2fl - Pi,o P2,3) [f-x]-'C ' (A4) 

on which that symmetry can be read off easily. 
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The 0[x3o] orthogonal matrix that effects the diagonahzation of [/-x]"^ can be used to 
re-define the integrations on fij and Q^. With this re-definition, the two Jacobians will 
compensate one another, so that keeping the same symbol for the re-defined f2's, under the 
integration over fij and Q^, one can proceed to the replacement 

ig^ni {[f-xYXo ^\^i9V^t — ^l (A5) 

where the ^^'s are the A^^— 1 eigenvalues of the matrix [/-Xso]) some of them, possibly zero. 

Now, relying on Theorem 3.2 in Ref. [3], and taking (1A3P and (lASP into account, it is 
possible to rewrite (lAip as 

1 /^"^ . ^ _ PAN Z^a = l?a oa (S"'' 06 

^ / d^i-.-da^Ee.) n le.-Ol ^ e-^^"^^"% (A6) 

■'^ 1 l<i<j<n V4l---4n 

where the delta-function accounts for the traceless property of any [f-x\ matrix, and where 
M is the normalization constant 

/+00 
dei---de„5(5^^.) n |e.-e,|e^^^-«^. (A7) 
'°° 1 i<«<i<" 

Of course, calculations can be continued further , but in (IA6I) it already appears that 
the possible occurrence of vanishing eigenvalues is not a problem, and that they should not 
contribute significantly. 



APPENDIX B: ESTIMATION IN SU(2) 

The quantities Y\l=i I-i ^3^5(1 — Xla (3^)^) ^^ed repeatedly in the text are simply solid 
angle factors, as can be seen immediately for = 3 of SU(2). There, one can make a variable 
change from 31, 32, 33 to A, 9, 0, where — 1 < A < +1, 31 = A sin 6* cos 0, 32 = A sin 6* sin 0, 
33 = A cos 6*. The Jacobian of the transformation d3id32d33 = JdAd^d0 is simple to obtain, 
J = sin 6, so that 

n r'd3^'^(l-$^(3'^)^) (Bl) 

/ + 1 pTV p27V 

d\ d9 d(j)5{l- X^)X^sm9 
-1 Jo Jo 

dO / d0sin^. 
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^0 



as expected. 

To illustrate how color fluctuations can reduce a coherent amplitude, consider the simplest 
SU(2) case of dUD, where 



J{n) = (271)-^ / d^t;e-^^-" ■ e™-^, (B2) 



Upon performing the angular integrations, and then integration over the magnitude of v, 
one obtains 



l_d_ ^ d_ 



(B3) 



The integral / d^f2 Ji(f2) e*"^ "^'^ then becomes 



^ n r^f^i^-iff) [l + ia{a\ + f)&^n^] e-^-^"^'''^', (B4) 

c=l "^"1 

and performing the final J d'^O J][(f2) on the result of (]B4p produces for the SU(2) form of 
(ITT]) the result 

d=i-^-i b 

■ {1 + laii [a^-G- (T^) + iai2 [a^ ■ G -l + i ■ G ■ a^) 
+zae3 (3 ■ G ■ 3) + (^«)'e4 ((7^ ■ G ■ 3) (3 • G ■ a^) 
+ M'e5 (3-G-3)' 



where ^1, ... ,^5 are numerical constants. The G""^ are gi-ven by that numerator function 
of ( l59l) . where the Mab depend upon the 3'^-components of Xso^ given by ( l58l) : and those 
3'^-components have been suppressed, for they require a separate, numerical integration. 
In SU(2), only the diagonal (Xg^ spin matrices can contribute to matrix elements between 
unchanged isotropic {i.e., color) states, whereas for SU(3) there would be two such matrices, 
Ag'^ and Ag^. 

Let us estimate the a 7^ effect by evaluating the '1' term of the curly bracket of ( IBSP ; 
and for this it is most convenient to consider an orthogonal transformation to diagonalize 
the real, ant i- symmetric G"*, by simultaneously transforming to a new set of variables 3'^, 
3'j,. Under such a transformation, the measures and form of the '1' terms contributing to 
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(IBSP are unchanged, but the exponential factor 3' • G" • 3' is simphfied because G' is diagonal. 
After converting to angular coordinates, let us simplify even further by suppressing the 0'-, 
(^'-dependence of that exponential, and merely calculate the 6*', 6' integrals, using 3' = cos 6' 
and 3' = cos 6\ 

I{a) = j ^ d^' J ^ dfe'-^"^', (B6) 

where a = g^{h)G'^^/r , and we assume that 6*33 7^ 0. These integrals are elementary and 
yield 

/(a) = - / dx . (B7) 

For large values of a, corresponding to small impact parameter, /(a) ~ 27r/a, and is damped 
in comparison to the corresponding integrals one finds for a ^ 0, /(O) = 4. This sort 
of damping is naturally to be expected when color fluctuations destroy the large impact 
parameter coherence, as stated in the text at the end of Section IIV[ 

APPENDIX C: GAUGE INDEPENDENCE 

The QCD Lagrangian can be expressed as 



■^QCD — -^gluon + -^quark + -^int (^1) 
1 

4 



where are gauge fields, F^^^ is the field strength with 

= d,Al - d^Al + gr'^AlAl, (C2) 
and A° are the color matrices of SU(3). Separate the gluon sector of Lagrangian into two 



parts as 



2l| 



-^gluon ^ iiv^ iiv 

\_ \ca. ca _|_ /-pia -pia ra ra \"| 

where f;j^ = d^^A^ - d^A^ as defined in Eq. Thus, £giuon = >Cg!uon + '^giuon with 

AO) _ _lfa va (CA\ 
gluon ^'-^lu'-^iu^ \^^> 



•^gluon ^\ fiu fiu tiv fiv J ' \^^/ 
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One can add and subtract a 'gauge-fixing' term to the gluon Lagrangian, which does not 
change its overall gauge invariance. 



(0) 

gluon 



1 1 

_fa fa Aa 



a\2 



•£ gluon l~' 



(0) 

gluon 



(C6) 
(C7) 



--(F;:,F;:,-f;j;j + -(5,A; 



a\2 
ji.) 1 



The QCD Lagrangian can then be written in terms of free -Cq^j-, and interacting 
parts as £qcd = -^qcd + '^qcd> where the free and interacting parts are 



(0) 

QCD 



quark 



(0) 

gluon 



c 



QCD 



-^int + -^1 



gluon 



(C8) 
(C9) 



1 

2C' 



respectively. 

The generating functional of QCD is 

1 



QCD 



i 6j' i 6f]' i Srj 



(CIO) 



where j^, rj^, and r;^ are gluon, quark and anti-quark sources, respectively. Following the 
conventional approach, either functional integral or Schwinger's Action principle 3|, the free 
generating functional with c'-^j.^ = C'-^il^^ + -Cquark is 



-2^o{j, V, V} = exp 



J 



(Cll) 



where the gauge field propagator is now defined by the gauge condition with a gauge pa- 
rameter ( as 



£(0) 
gluon 



+ 



(C12) 



a jrab 



A. 



\ / fjLV 
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and 



ab 



-6 



ab 



1 



- 1 d.d. 



After rearrangement, one finds 



J J '-'c J . g 2 J «A 



(C13) 



(C14) 



■e 



i _2_ ^ _£_ ■ r - a 
g 2 J J J . e 2 J 5A SA 



gL[A] 



(S) 



where A^{x) = J dy Dc^u {x — y) j^{y)- The exponential factor involving >Cgiuon can be cast 
into the form 



ifC, [A] 

p J gluoilL J 

exp j (f^^f;:, - f;j;j + ^ | (s.A^D^j 



(C15) 



where Eq. (fT7|) is used. The Xni.-independent factor becomes 



(C16) 



2 , ^l^''' 



1 



The generating function of QCD becomes 

= e'^ ■' ■ e 2 J SA SA 



(C17) 



■e ' 2 



,,1-1 pL[A] 



(S) 



Except for the expansion of the closed-fermion-functional L[A], the gauge field dependence 
in the exponent is at most quadratic; however, an expansion in powers of using a 

modified Fradkin representation for each L[A] and Gc[A], generates a totally quadratic 
A-dependence. 
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For the QQ or QQ scattering, one will encounter 



e 2 J SA '-^t: ■ SA ■ 



Gl[A] Gl[A] e^[^] e+t/x-[f+9M^]+f /A-D, 



(0-1 



(C18) 



A-*0 



Under the eikonal and quenched approximations, the coefficients of linear and quadratic 
y4^-dependent terms are (c/. Eq. fl25|) ) 



and 



ab 



(C19) 



(C20) 



respectively, and where T^fj, and TZ^^ come from the eikonal approximation of the Green's 
function of the quarks or anti-quarks. The linkage operation can be worked out as 



e-2/ A-dF-A . e+l/A-'c-A+i/A.Q 







i:>P-( 1-/C-D 



■Q 



(C21) 



The kernel in the quadratic term of is 



(1 - /C ■ Bif 



{gf -x)^ ■ 



(C22) 



The result is independent of the gluon (gauge field) propagator. The derivation is valid for 
arbitrary relativistic gauge conditions. 

APPENDIX D: LIST OF ABBREVIATIONS 



The following abbreviations have been used freely throughout the text. 
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CM 


Center of Mass 


ETCRs 


Equal-time Commutation Relations 


FI 


Functional Integral 


GF 


Generating Functional 


GI 


Gauge- Invariant 


LC 


Lorentz Covariant 


MGI 


Manifestly Gauge Invariant 


MLC 


Manifestly Lorentz Covariant 


NVM 


Neutral Vector Meson 


OE 


Ordered Exponential 


Q 


Quark 


Q 


Anti-quark 


QA 


Quasi-Abelian 


QFT 


Quantum Field Theory 


RHS 


Right Hand Side 
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